Novel interface-selected waves and their influences on wave competitions 



Xiaohua Cui 1 , Xiaoqing Huang 1 , Zhoujian Cao 2 , Hong Zhang 3 , and Gang Hu^ 
Department of physics, Beijing Normal University, Beijing 100875, China 
2 Institute of Applied Mathematics, Academy of Mathematics and System Science, 
Chinese Academy of Sciences, 55, Zhongguancun Donglu, Beijing 100080, China and 
3 Zhejiang Institute of Modern Physics and Department of Physics, Zhejiang 
University, Hangzhou 310027, China 
(Dated: June 20, 2008) 



The topic of interface effects in wave propagation has attracted great attention due to their the- 
oretical significance and practical importance. In this paper we study nonlinear oscillatory systems 
consisting of two media separated by an interface, and find a novel phenomenon: interface can select 
a type of waves (ISWs). Under certain well defined parameter condition, these waves propagate in 
two different media with same frequency and same wave number; the interface of two media is trans- 
parent to these waves. The frequency and wave number of these interface-selected waves (ISWs) are 
predicted explicitly. Varying parameters from this parameter set, the wave numbers of two domains 
become different, and the difference increases from zero continuously as the distance between the 
given parameters and this parameter set increases from zero. It is found that ISWs can play crucial 
roles in practical problems of wave competitions, e.g., ISWs can suppress spirals and antispirals. 

PACS numbers: 



The behaviors of waves around interface between two media have attracted continual and great interest [1-7]. In 
linear optics we are familiar with the problems of wave reflection and refraction, which are predicted analytically. 
However, for nonlinear systems the interface-related behaviors become much more complex and diverse, and much 
less known. 

The problems of wave propagation in linear and nonlinear media have also attracted considerable attention. Re- 
cently, new observations of inwardly propagating waves have stimulated considerable interest in this field. For several 
centuries, scientists have known waves propagating forwardly from wave source only, called normal waves (NW). In 
recent years, different types of waves propagating toward the wave source (called here antiwaves, AW) have been 
observed in both linear optics [1, 8] and nonlinear oscillatory systems [9-14]. These new phenomena introduce com- 
pletely new topics of the interface problem. For instance, novel phenomenon of negative refraction has been reported 
in both linear optics [1, 8] and nonlinear oscillatory systems [3,7]. 

In the present paper, we find another completely new nonlinear interface-phenomenon: interface of two different 
media can generate waves, called here interface selected waves (ISWs). In a well defined parameter surface the 
frequency and wave number (also wave length) of ISWs are identical in two media with different parameters, and they 
can be predicted analytically. Away but near this surface ISWs still exist, though the above analytical predictions 
are no longer available and wave numbers of ISWs in the two domains are no longer identical. By varying parameter 
away from this surface continuously, wave numbers change continuously, and the difference of wave numbers in two 
domains also increases from zero continuously. It is found that ISWs play crucial roles in practical problems of wave 
competitions in oscillatory systems, e.g., in suppressing spirals and antispirals. 

We consider the following bidomain reaction-diffusion system 
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where D(^i) isamxm matrix with constant elements depending on 7.;. The function f and diffusion matrix D in the 
two domains are identical because the same reaction-diffusion processes occur in both sides. On the other hand, the 
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dynamical evolutions in each side may be different due to different control parameters (/xi, v\, 71) and (/X2, v 2 , 72)- 
We represent the interface between the two domains by /, then the following boundary conditions are required on I 



Ui{/}=U2{/} ' d^r } = d^r } (2) 

where indicates a derivative of U, over space variable along the direction perpendicular to the interface /. We 

assume that Ui = U2 = is a stable point of Eq.(l) and fii(vi) controls (not controls) the linear terms of Uj for the 
reaction parts. Hopf bifurcation with frequency lo is supposed to occur at parameters /xi = /x 2 = /x . Moreover, we 
assume further that both /xi and /x 2 are slightly beyond the Hopf bifurcation point 



^( Mi W- M o (j ' ) ) 2 «l, * = 1,2 (3) 

At fii, XJi performs periodic oscillation of frequency Uj, and u>i approaches to luq as fii reduces to /xo (i = 1, 2). Under 
condition(3) Eq.(l) can be reduced to amplitude equations, i.e., the bidomain complex Ginzburg-Landau equations 
(BCGLE) by the approach standard for the derivative of single-domain CGLE [15, 16]. 

-^-= ai (l- itl 1 )A 1 - 6i(l + ia 1 )\A 1 \ 2 A 1 + ci(l + */?i)V 2 Ai (4a) 



dt 



= a 2 (l - ^fi 2 )A 2 - 6 2 (1 + *a 2 )|A 2 f A 2 + c 2 (l + i[3 2 )V z A 2 (4b) 

Mi dA 2 

A 1{1} = A 2{1} , - = 4c 

where (a^, fii), Qj) and (cj, are related to /i., — /i.o, (/Ltj — /Lio, I'i) and (/Xj — /Lio, 7,), respectively. The continuity 
conditions Eq.(4c) can be derived from condition (2) because the transformation from f/i to Ai is exactly the same as 
that from U 2 to A 2 (in the two domains, amplitude equations are derived at a common Hopf bifurcation point with 
an identical linear matrix at /xo), and the transformations from Ui to Ai (i — 1, 2) are determined only by this linear 
matrix. In Eq.(4) A\ and A 2 are complex variables in each side of the interface. With scaling transformations we can 
fix ai, 61, Ci and fii, and the remaining 8 parameters are irreducible for BCGLE systems. In the following study we 
will set ai = 61 = Ci = 1, fii = for numerical simulations without mentioning and all the theoretical formulas are 
given generally for 12 parameters. Without the interface interaction, the two media have their single-domain planar 
wave solutions [2, 17] 



M,r. / 1 = J(±(oi - Cifc?))^*"-*"'), < fc 2 < I (5a) 



u>i = a,i(ai + fii) + a(0i - ai)k 2 , a t , b t , c, L > (5b) 
Waves in media Mi and M 2 are classified to NWs and AWs under the conditions [3, 12, 14] 

NWs : cj,a,(fij + a{) < ; or, w,aj(fij + ai) > and \u>i\ > |aj(fii + cei)\ (6a) 



AWs : ujiai(rii + c^) > and < |aj(fij + CKj)| (6b) 

By AWs we mean waves with negative phase velocity, while both NWs and AWs have positive group velocities [11, 
12, 14]. Now we focus on the interface related problems, and start from a one-dimensional (ID) BCGLE system. We 
are interested in the problem how the interface can significantly influence the system dynamics. In Fig.l(a)-(c) we 
study the system evolutions at two different parameter sets with random initial conditions, and find characteristically 
different features in the asymptotic states. The most significant and new observation is given in Figs. 1(a) and 1(b) 
where we find homogeneous planar waves moving in both media from right to left with a constant velocity and 
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transparently crossing the interface. These homogeneous running waves originate from the interface (see Fig. 1(b)), 
therefore, called interface selected waves (ISWs). The phenomenon of Fig. 1(a) is surprising. With two media having 
different control parameters we intuitively expect that the waves in two media must have different wave numbers (even 
if both sides may have the same frequency). This common feature is clearly seen in Fig. 1(c) where we observe uniform 
bulk oscillation in the right medium and waves propagating from left to right in the left medium, clearly manifesting 
the interface /. However, in Fig. 1(a) waves propagate seemly in a homogeneous medium without feeling any difference 
of Mi and M 2 ; the interface is transparent to the waves. Moreover, the realization of these running waves is stable 
against different initial perturbations. This characteristic phenomenon can never exist in linear systems, and has 
never been observed so far in nonlinear systems. 

It is interesting that we can predict the frequency and wave number of ISWs explicitly and exactly under certain 
parameter conditions. For case of Fig. 1(a) we can determine the frequency and wave number by the following simple 
requirements: 

wi(fci) = w 2 (fe 2 ), fci = k 2 , | Ai| = \A 2 \ (7) 
Inserting Eqs.(7) into Eq.(5) we obtain a unique set of solutions ojj, kj 

k 2 = k 2 = k 2 = ai(tti + »i)-a 2 (tt2 + » 2 ) (8a) 
1 2 1 c 2 (/3 2 -a 2 )-ci(/3i-ai) 

ai(ai + Oi)c 2 (/3 2 - a 2 ) - Ci(/3i - a 1 )a 2 (a 2 + fi 2 ) . 

lu 1 = uj 2 = u>j = — — (8b) 

c 2 (p 2 - a 2 ) - ci(/?i - ai) 

And these solutions are exact in the parameter surface 

for cib 2 - c 2 bi j= : if and only if — — — = (8c) 

ci& 2 -c 2 &i c 2 {p 2 - a 2 ) - ci(j»i - ax) 

for c\b 2 — c 2 bi = : if and only if a\b 2 — a 2 b\ = 

which is obtained by inserting Eq.(8a) into Eq.(5a), and by the identifying \A\\ = \A 2 \. It can be easily confirmed 
that on the parameter surface Eq.(8c) the planar wave solution Eq.(5) with frequency and wave numbers given by 
Eqs.(8a) and (8b) are exact solution of BCGLE. Moreover, the predictions of Eqs.(7) and (8) agree exactly with the 
numerical results of Fig. 1(a). 

In Figs. 2(a) and 2(b) we specify the surface of Eq.(8c) in some parameter planes where solid lines represent the 
parameters satisfying condition (8c). In Figs. 2(c) and 2(d) we vary parameters along the solid line of Fig. 2(a) and 
numerically compute ID BCGLE. we compare the numerical results (empty cycles and triangles) with theoretical 
predictions of Eqs.(8a) and (8c) (solid lines), and find: (i) ISWs exist in a large area of surface Eq.(8c); (ii) the 
predictions of Eqs.(8a) and (8b) coincide with numerical results exactly (within computation precision). In Fig. 2(e) we 
fix a parameter set on the surface (black disk T) and present the asymptotic pattern evolution of the BCGLE system. 
It is clearly shown that ISWs, with the wave number and frequency predicted by Eqs.(8a) and (8b), asymptotically 
control the entire bidomain during their propagation. For demonstrating the possibility of observation of ISWs in 
experiments, we study a reaction-diffusion model: bidomain Brusselator. In Fig. 2(f), we show ISWs of this chemical 
model, satisfying all conditions of Eqs.(7). 

The solutions of Eqs.(8a) and (8b) are exact for BCGLE only on the parameter surface of Eq.(8c). Slightly away 
from this surface Eqs.(8a) and (8b) can no longer predict the wave numbers and the frequencies of ISWs exactly. In 
Figs. 3(a) and 3(b) we compare the theoretical predictions of Eqs.(8a) and (8b) with numerical results of frequencies 
and wave numbers by varying parameters along the dashed line in Fig. 2(a). We find: (i) The solutions of Eqs.(8a) 
and (8b) are not exact; the wave numbers in the two domains deviate from each other (about 7.7% difference in 
Fig. 3(d)); (ii) However, slightly away from the surface Eq.(8c), the feature that the interface generates waves is still 
clearly observed (compare Fig. 1(b) with Fig. 3(c)), i.e., ISWs still exist; (iii) By continuously increasing the parameter 
distance from condition Eq.(8c), the deviation of the numerical results from the theoretical predictions Eqs.(8a) and 
(8b) increases continuously from zero too, and for small parameter deviation the solutions Eqs.(8a) and (8b) can still 
be used for predicting frequency and wave numbers with very good approximation. In Figs. 3(c) and 3(d) we show 
ISWs for a parameter set away from the surface Eq.(8c) (Disk Q in Fig. 2 (a)). It is clear that even away from surface 
Eq.(8c) the waves of Fig. 3(c) are generated by the interface in the similar way as in Fig. 1(b) though we have k\ = k 2 
in Fig. 1(b) but k\ ^ k 2 in Fig. 3(c). Thus the waves in Fig. 1(a), 1(b) and Figs. 3(c), 3(d), have obviously the same 
interface-selected nature which is essentially different from the waves of Fig. 1(c). Similar ISWs with k\ ^ k 2 can also 
be observed in bidomain Brusselator. In Figs. 3(e) and 3(f) we take parameter set far away from that in Fig. 2(f), and 

can still observe ISWs. Here the wave numbers in the two sides have slight difference ( ^^^ ~ 6.77%). 
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There are some necessary conditions for ISWs to appear. Let us analytically specify some of these conditions under 
Eq.(8c) (Eqs.(8a) and (8b) are exact solutions of uj and kj). From Eqs.(8a) and (5a) we have an obvious necessary 
existence condition for ISWs, i.e., 



■ 2 ai(ai + fli) - a 2 (a 2 + Q 2 ) . a { 



0<kf = -^ ^ ^ y<-, * = 1,2 (9a) 

ci\pi - a 2 ) - ci(/?i - an) a 

which should be satisfied because wave number ki must be real. If this condition is violated, there is no physically 
meaningful solutions of ki and u>i, and thus no ISWs can be observed. This is the case of Fig. 1(c). ISWs are 
generated by the interface and propagate along a certain direction. Therefore, the waves must forwardly propagate 
in one domain. Precisely, ISWs are NWs in the left (or right) domain while AWs in the right (left) domain for waves 
propagating from right (left) to left (right), and this requires another parameter condition 

AWs : wjaj(fli + a») > and < |a,(Qj + aj)| (9b) 
NWs : u)iaj(Clj + cq) < or u)iaj(Qj + cq) > and > \ai(Qj + cq)\ 
i = 1, i = 2 or i = 2, i = 1 

In order to provide an idea how these conditions influence the existence of ISWs, we show Fig. 4, where one can 
numerically observe ISWs in the regions enclosed by disks, called "ISW" regions. In Fig. 4 dashed-dotted lines with 
kj = and a, = i — 1 or 2 are the boundaries of "ISW" theoretically predicted by Eqs.(9a)and (9b), respectively. 
In "No ISW" regions ISWs do not exist due to violations of conditions of Eq.(9a) or (9b). In the regions "Unstable" , 
ISWs exist while waves with wave number ki arc unstable due to Eckhaus instability, and there ISWs cannot be 
numerically observed. Figures 4(a) and 4(b) are plotted in a small parameter surface under the condition of Eq.(8c). 
Similar structure of distributions "ISW", "No ISW" and "Unstable" regions can be observed when parameters are 
varied slightly away from the set of Eq. (8c) . 

Though the above investigations are made for ID bidomain systems, the observations of ISWs exist generally for 
high-dimensional systems. In 2D oscillatory systems, much richer types of waves, including spirals and antispirals, 
can be self-sustained, and wave competitions become an important issue. Now we explore how ISWs play crucial roles 
in wave competitions. We consider a 2D BCGLE system with an interface line IV in between. Without the interface 
interaction, Mi supports normal spirals (Figs. 5(a), (d), (g)) and M 2 supports antispirals ( Figs. 5(b), (e), (h)). With 
the interface interaction we find characteristically different results of wave competitions. Fig. 5(c): the antispiral wins 
the competition and dominates the system with frequency w 2 ; Fig. 5(f): the spiral wins; Fig.5(i): ISWs win and 
dominate the two domains. The reasons why we can observe so diverse results in similar competitions between spiral 
and antispiral waves can be completely understood based on the analysis of ISWs. 

For explaining the results of Fig. 5 we briefly introduce some known conclusions on wave competitions in oscillatory 
systems. If competitions occur in a homogeneous medium, the results are [3, 18]: 

NWs against NWs : faster waves win (10a) 



AWs against AWs : slower waves win (10b) 

NWs against AWs : NWs win (10c) 

With the competition rules Eq.(10) and the analytical results of Eqs.(7)-(9) we can fully understand and predict the 
diverse results of Fig. 5. 

Inserting the parameters of Fig. 5(c) into Eq.(8) we have u>i = 0.636. Considering conditions Eq.(6) we conclude 
that ISWs are NWs in Mi and AWs in M 2 (note, the interface is the source of ISWs). The frequencies of the spiral in 
Mi and the antispiral in M 2 are lo\ = 0.288 and oj 2 = 0.626, respectively. According to conclusion (10) ISWs win the 
competition in Mi against the spiral while losing the battle in M 2 against the AW spiral. Therefore, the antispiral 
waves of frequency w 2 finally dominate the whole system. The parameters of Fig. 5(f) do not satisfy condition Eq.(9a), 
and no ISWs can be generated. In Figs. 5(d) and (e) we observe lo\ — —0.0567 and w 2 = 0.151. According to Eq.(6) 
waves of o>i(w 2 ) are NWs (AWs) in both Mi and M 2 . On the basis of (10c), the spiral of frequency u\ wins the 
competition. The most interesting observation is given in Fig.5(i) where we have u>i — 0.932. The frequency of the 
spiral(antispiral) in Mi (M 2 ) is lo\ — 0.293 (a> 2 = 1.037). Therefore, ISWs win both competitions against the spiral in 
Mi (condition (10a)) and against the antispiral in M 2 (condition (10b)). The asymptotic state is ISWs in 2D system 
where ISWs suppress both spiral and antispiral in Fig.5(i). 

In conclusion, we investigated the role played by interfaces. A new type of waves, interface selected waves (ISWs) 
were found in bidomain systems with one medium supports AWs and the other NWs. When control parameters are 
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on a well denned parameter surface ISWs propagate with analytically predictable same frequency and wave number 
in two media with different parameters. When the parameters are away but near this surface, ISWs can be also 
observed, of which the frequency and wave numbers can be located approximately. These waves are selected by 
interfaces between two media, and some necessary conditions for observing these ISWs are specified. These ISWs play 
important roles in wave competitions. For instance, under certain conditions ISWs can suppress spiral and antispiral 
waves in both media. These roles are important in practical applications. Experimental realizations of ISWs in 
chemical reaction-diffusion systems are strongly suggested, based on the well-behaved ISWs of Fig. 2(f), Fig. 3(e) and 
Fig. 3(f) computed for a chemical reaction-diffusion model. 
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FIG. 1: (a)-(c) Contour patterns of HeAi of a ID BCGLE system with interface /. Domains Mi, (i = 1, 2) have length L and 
parameters ai, f3i, and a,i = bi = Cj = 1, Q.i =G, i = 1, 2. Numerical simulations are made with space step dx = 0.5, time step 
dt = 0.005, and L — 200. No-flux boundary condition, randomly chosen initial conditions and the above time and space steps 
are used in all figures for numerical simulations unless specified otherwise, (a) ot\ — 0.1, (3\ = —1.4, Q2 = —0.2, 02 = 1.4. 
Interface-selected waves (ISWs) homogeneous in Mi and M2 are observed, and the interface is transparent to ISWs. (b) The 
same as (a) with early time evolution plotted. It is clearly shown that ISWs originate from the interface, (c) qi = 0.1, /3i = 
— 1.4, o?2 = 0.5, P2 = 1-4. Mi and M2 support different regular waves, clearly manifesting interface /. 




FIG. 2: (a) Surface Eq.(8c) (Solid line) is plotted in 0,2 — 62 plane with Q2 = 0, C2 = 2.0, Qi = 0.6, /3i = —1.4, Q2 = 
0.1, /?2 = 1.2. (b) The same as (a) with surface (8c) plotted in 012 — (12 plane, 62 = 1.778. Point E corresponds to the 
boundary k 2 = 0. (c) Frequency of ISWs which are selected along the solid line of (a). Theoretical prediction Eq.(8a) 
(solid line) coincides with numerical simulation for ID BCGLE (empty circles and triangles) perfectly, (d) The same as 
(c) with wave numbers ki = &2 plotted. Agreement between theoretical prediction Eq.(8a) and numerical results is also 
confirmed, (e) ISW pattern obtained by using the parameter set a2 = 1.8, 62 — 1.778 (point T in (a)), (f) The same 
as (e) with contour pattern of v variable of Brusselator RD which is numerically computed for ID chain. The system is: 
= ai-ibi+l+^u+tl+ai^v+Sut^u, ^f- = biU-(l+<Ti)u 2 v+6 vi V 2 v, i = 1,2, ai = 1.0, 61 = 2.24, 71 = o\ = 0.0, S ul = 
2.31, S v i = 2.17; a 2 = 1.02, b 2 = 2.2624, 72 = 0.02, cr 2 = 0.01, S u2 = 0.95, S v2 = 2.47; dx = dy = 0.5, dt = 0.0025, L = 300. 
ISWs with identical w and k are observed. 
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FIG. 3: (a) (b) The same as Figs. 2(c) and (d), respectively, with parameters varied along the dashed line of Fig. 2(a). Numerical 
simulations are made for ID BCGLE. Now deviations between theoretical results of Eqs.(8a), (8b) and the numerical results are 
observed. Deviation increases as parameters vary away from the surface Eq.(8c). (c) (d) The same as Fig. 2(e) with different time 
intervals plotted, respectively, in which parameters are taken away from the solid line of Fig. 2(a) (point Q, a 2 = 1.6, b 2 = 1.778). 
Now ISWs are still observed while wave numbers in the two sides are slightly different jH^iy w 7.71%, Ak = \ki - k 2 \). (c) 
(f) The same as (c) and (d), respectively, with ID Brusselator chain computed. The parameter set is considerable different 
from that of Fig.2(f): 01 = 1.0, 61 = 3.2, 71 = ai = 0.0, <5„i = 1.0, 8 vl = 0.5; a 2 = 1.1, b 2 = 3.2, 72 = 0.1, a 2 = 0.0, 8 u2 = 
1.0, S v2 — 2.5. Now the wave numbers of the two sides are also slightly different ~ 6 - 77 %)- 
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FIG. 4: The distributions of different types of waves in (ai, a 2 ) parameter planes for different sets of (/3i, f3 2 ). Black disks 
represent the boundaries of ISW regions ("ISW" regions) identified by direct numerical simulations of ID BCGLE. In "No ISW" 
regions, ISWs do not exist due to the violations of condition Eq.(9a) (boundary kj = 0, i.e., oti = a 2 ) or condition Eq.(9b) 
(boundary a, = /3i, i = 1 or 2) (both presented in dashed-dotted lines). In the region "Unstable", both conditions Eqs.(9a) and 
(9b) are satisfied, but waves with the given ki are unstable due to the Eckhaus instability, ai = 6» = Cj = 1, fi» = 0, i = 1, 2. 
(a) /3i = — 1.4, /3 2 = 1.4, (b) f3i = 2.5, /3 2 = 0.5. Black triangles A and B in Fig. 4(a) represent the parameter sets used in 
Figs. 1(a), (c), respectively. 
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FIG. 5: Wave competitions between spiral, antispiral and ISWs in 2D BCGLE with interface II' . at = bi = Cj = 1, Qi = 0, i = 
1,2. (a)(d)(g) Spirals in Mi medium. (b)(e)(h) Antispirals in M2 medium. Snapshots in (a)(b), (d)(e) and (g)(h) are used as 
the initial conditions for the dynamic evolutions of (c), (f) and (i), respectively. (c)(f)(i) The asymptotic states of bidomain 
systems with interface II . (a)(b)(c) Qi = 0.2, f3\ = 2.0, Q2 = 1.0, 02 = —0.5. Antispiral initially in Ah dominates the system 
in (c). (d)(e)(f) ai = 0.2, /3i = -2.0, a 2 = 0.5, (3 2 = —1.4. Spiral initially in Mi dominates the system in (f). (g)(h)(i) 
ai = 0.1, Pi = 3.2, o?2 = 1.2, P2 = 0.0. In (i) ISWs suppress both spiral in Mi and antispiral in M2, and dominate the whole 
system. 



